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Abstract 


We derive the wave equation obeyed by electromagnetic helds in curved space- 
time. We find that there are Riemann and Ricci curvature coupling terms to the 
photon polarisation which result in a polarisation dependent deviation of the photon 
trajectories from null geodesics. Photons are found to have an effective mass in an 
external gravitational field and their velocity in an inertial frame is in general less 
than c. The effective photon mass in the Schwarzschild metric is = (2GM 
and near the horizon it is larger than the Hawking temperature of the blackhole. 
Our result implies that Hawking radiation of photons would not take place. We also 
conclude that there is no superluminal photon velocity in higher derivative gravity 
theories (arising from QED radiative corrections), as has been claimed in literature. 
We show that these erroneous claims are due to the neglect of the Riemann and Ricci 
coupling terms which exist in Einstein’s gravity. 



A standard result of Einstein’s gravity is that the trajectories of all massless 
particles are null geodesics. A question worth examining is whether there is a devi¬ 
ation from the null geodesics if the particles have a spin i.e. due to the interaction 
of spin with the Riemann and Ricci curvatures of the gravitational helds. In this 
letter we have studied this question for the case of photon propagation in a curved 
background. Starting with the action ^/—gF^iy for electromagnetic helds in a 
gravitational held, we derive the wave equation for electromagnetic held tensor F^y, 
which turns out to be of the form (Eddington [1] , Noonan [2]) 

- R\F,p + = 0 (1) 


We see that the photon propagation depends upon the coupling between the Riemann 
and Ricci curvatures and the photon polarisation. This leads to a deviation of the 
photon trajectories from the null geodesic by amounts proportional to the Riemann 
and Ricci curvatures. The photon trajectories in the geometrical optics limit are 
described by the following generalisation of the geodesic equation 


(RX^ dX^ dX^ 

ds^ ds ds 


:V“ [R^\Jp, + R\U 


RRfxp] X 


P 


( 2 ) 


The nonzero right hand side of the modihed geodesic equation (2) leads to a polari¬ 
sation dependence in the gravitational red shift, bending, and Shapiro time delay of 
light. The redshift in a photon’s wavelength propagating in the gravitational held of 
a spherical mass (or a blackhole) is given by 

A2 /l-2GM/r2y/' 

Ai “ l^l-2GM/rJ 

The extra factor in the curly brackets is the correction to the standard red shift 
formula due to Riemann curvature coupling. This deviation is measurable when 
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the photon wavelength is of the same order as square root of the magnitude of the 
Riemann curvature. Near the horizon of blackholes, where the curvature can be 
large, this deviation is signihcant. We derive the dispersion equation of a photon in 
the Schwarzschild metric from the wave equation (1) and hnd that in a local inertial 
frame the dispersion relation is given by 

2 T2 ‘^GM 
ui - k = ^3 

This implies that in the gravitational held of a blackhole, photons acquire an effective 
mass 

/2GM\^G 

) 

At the horizon r = 2GM, the effective photon mass is = 1/{2GM). Since the 
photon mass is larger than the Hawking temperature Th = I/SttGM, our result im¬ 
plies that Hawking radiation of photons would not take place. 


A curious phenomena discovered by Drummond and Hathrell [3] is that in 
higher derivative gravity which arises by QED radiative corrections, the photon ve¬ 
locity in local inertial frame can exceed the velocity of light in the Minkowski space. 
Due to the coupling of the electromagnetic helds with the Riemann and Ricci ten¬ 
sors in the Lagrangian, it was claimed that the photon velocity in the Schwarzschild, 
Robertson- Walker, gravitational wave and deSitter backgrounds is larger than the 
flat space velocity c. This result was extended by Daniels and Shore to charged [4] 
and rotating blackholes [5] with the same conclusions. Latorre et al [6] have shown a 
universal relation between the velocity shift of photons to the energy density which 
generates the background metric, and Shore [7] has related the velocity shift to the 
coefficients of conformal anomaly. Lafrance and Myers [8] interpret these results as 
the breakdown of the Equivalence principle in higher derivative gravity and Dolgov 
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and Khriptovich [9,10] derive this result from field theoretic dispersion relations. Fi¬ 
nally Mende [11] has proposed this effect as a test for string theories of quantum 
gravity. 


In this letter we show that claims of superluminal photon velocity are due to 
the neglect of the Riemann coupling terms in the wave equation which arises from the 
minimal Lagrangian. We hnd that in Einstein’s gravity the photon velocity 

in a Schwarzschild blackhole metric is given by 

c 


and for the Friedman-Robertson-Walker metric it is given by 

c 

+ ISj - 3p)) ' 


(7) 


Therefore the velocity of photons in a gravitational held is less than c in Einstein’s 
gravity (as long as the strong enegy condition p > 3p is satished. We derive the wave 
equation by including higher derivative terms and show that Riemann coupling terms 
from the higher derivative gravity are always smaller in magnitude than the Riemann 
term that already exists in Einstein’s gravity. This analysis shows that the photon 
velocity never exceeds c, whether in Einstein’s theory or by the inclusion of higher 
derivative terms in the Lagrangian. 


The interaction of electromagnetic helds with gravity is given by the action 

F = y d^x^gF^^F^- (8) 

(We use the cfonvention c = 1, signature —2 and Greek letters denote spacetime 
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indices 0 — 3). From (8) we obtain the equations of motion 


= 0 (9) 

Equation (9) with the Bianchi identity 

+ Vi/Fa/^ + VxFfj_i, = 0 (10) 

gives the Maxwell’s equation in curved background. Operating on equation (9) by 
the Va and using the Bianchi identity (10) the wave equation may be obtained as 

V^V^FjyA + Va (V^Fjy^) + [V^, V;/] Fa^ 

( 11 ) 

- [Va,V^]F^, = 0 

The second term vanishes owing to (9). Using the identity for the commutator of 
covariant derivatives: 


[V^ V,] F,^ = + RpxF^ 


and the circular identity 

Rp\u^ T RpupX T RppXv 0 
the wave equation (11) reduces to the form 


( 12 ) 


( 13 ) 


+ R\F,, - R\F>„ = 0 ( 14 ) 

The Riemann and Ricci curvature coupling terms to the photon polarisation (or spin) 
give rise to the polarisation dependent deviation of photon orbits from null geodesics. 


Photon trajectories are described by the eikonal solutions of the wave equation 
(14) of the form 


F = 

liv 


— JS{x) 


f, 




( 15 ) 
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where the phase S varies more rapidly in spacetime than the amplitude f^,y. The 
wavenumber of the photon trajectories is given by the gradient of the phase 


= V,S 


In the geometrical optics approximation where 


(16) 


the wave equation (14) can be written in the form 

- .Jp, - R\fup + R^Jxp = 0 


(17) 


(18) 


The dispersion relation may be written as 

A'2 = (a» (19) 

Operating on (19) by V", the L.H.S. is 


V“/s:2 = (20) 

where we have used the identity which follows from the dehnition of 

as a gradient. Light rays are dehned as the integral curves of the wave vector 
i.e. the curves x^{s) for which = K^. Substituting = dX^/ds in (20) we have 
for the L.H.S. 


V«A- = 2^V, 


ds 


ds 


= 2 




ds^ 


ds ds 


( 21 ) 


Using (19), (20) and (21) we obtain the modihed geodesic equation for photon tra¬ 
jectories 


*2 ds ds 2^ ■'A /«. R xU + R Ms) X 


( 22 ) 
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To obtain the photon velocity in the curved space, we can use the dispersion relation 
(18) directly. Of the six components of in equation (18), only three are indepen¬ 
dent owing to the Bianchi identity (10). Choosing the components of the electric held 
vector Ei = fd as the independent components we have from the Bianchi identity 
( 10 ) 

Kofij + Kifjo + Kjfoi = 0 (23) 

Using (23) to substitute for in terms of the electric held components fjo and foi 
in the wave equation (18) we obtain 

Y+A'^/O + R\ + R' of) o.f + R’. 

which is the wave equation obeyed by the three components of the electric held vector 
(in (24) Latin indices are over three spacelike components and repeated indices are 
summed over). To simplify the notation we write (20) as 

{KHl + e{)foj=0 (25) 

where 

€l=[R\ + R\!Fyi+(^-2R‘’\, + iR“ (26) 

Of the three components of the electric held vector in the electromagnetic waves, 
only the transverse waves are observable in an asymptotically hat space. To obtain 
the equation for the transverse helds we hrst diagonalise the matrix el to give the 
equation for the three normal modes: 

(A:' + Q)/oi = 0 (y = 1,2,3) (27) 

where e* are the eigenvalues of the e* matrix. The three equations (27) can be sepa¬ 
rated into the equation for the transverse helds 

/y = - n‘n) /„ (28) 
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and the longitudinal fields 


foj^ = n^njfoi (29) 

where n* = are the components of the unit vector along the direction of propaga¬ 
tion. The wave equation for the transverse photon is 

+ ei (Sj - foj = 0 (30) 

and for the longitudinal photon the wave equation is 

+ e*) foj = 0 (31) 

Non-Rotating Black Holes (Schwarzschild Metric): 




For tangential trajectories, n = (0,1, 0) and the equations for the transverse helds Ei 


and yield: 


( K^ + AGM/r^ 
[ 0 


0 

2GM/r 



The dispersion relation is (^K‘^ + {k"^ — = 0 

to the propagating mode is 


I ) = 0 (37) 

and the root corresponding 


=oj^-k 
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(38) 


and the velocity of wave propagation is given by 

du 1 


T) _ 


' = 


dko 


(l + w (^)) 


1/2 


< 1 


(39) 


Thus for both radial and tangential trajectories the photon velocity is subluminal and 
the effective photon mass is given by (36). From the dispersion relations (34) and 
(38) for the radial and tangential photons we see that a photon with frequency u (in 
the inertial frame) reaches a terminal radius = {2GM at which point its 
wavenumber becomes zero. A hnite energy photon falling into a blackhole therefore 
would never encounter the singularity at the origin of the blackhole. 


At the horizon of the blackhole the effective photon mass is = ^GM. Com¬ 
paring this with the Hawking temperature of the blackhole = ^ttGM, we hnd that 
the photon mass is larger than the Hawking temperature. This implies that Hawking 
radiation from blackholes cannot be in the form of photons. 


In the coordinate frame the dispersion relation obeyed by = (u, Kkg, kp) 
obtained from the wave equation (30), for a radial trajectory is given by 


1 - 


2GM\ 


-1 
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(40) 
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Since cu is a constant of motion (the metric being stationary), the wavenumber kr 
changes with r. The redshift of the wavelength A = {grr)^^'^ /K is given by the 
expression 

A(r2) _ /l-2GM/r2y/' 

A(ri) l^l-2GM/rJ 

The extra factor in the curly brackets is the correction to the standard redshift for¬ 
mula due to the coupling of the Riemann curvature to the wave vector. The correction 
is signihcant when the photon wavelength is comparable to the Riemann curvature. 
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Expanding Universe (Friedmann-Robertson-Walker Metric): 


For wave propagation in a general homogeneous, isotropic universe described 
by thge line element ds'^ = +dt^ —R^{t) {dr"^ -|- r'^dO'^ + sin^ the non-zero com¬ 

ponents of the e* matrix (26) are 


e 


1 

1 



-2 


'R R^' 


SttG 


(P - 3p) 


(42) 


where R is the scale factor, and p and p are the pressure and energy density respec¬ 
tively. 


The dispersion relation for the transverse photons is 

SttG 


cu" - ki = 


(P - 3p) 


(43) 


and the photon velocity is given by 



Therefore we find that the photons are not superluminal as long as the strong energy 
conditon p > 3p is satished. In the radiation dominated era when p = 3p the photon 
velocity is 1. One may turn this argunment around and use the axiom from special 
relativity that photons of any polarisation cannot exceed the flat space velocity c when 
measured in an inertial frame, to show that the strong energy condition p — 3p > 0 
must be obeyed for any form of matter. Consider a space time region hlled with some 
matter with equation of state p = qp. The e* matrix in that region of spacetime 
is given by e* = —Sj{8nGp/3) (1 — {3/q)) (assuming the matter distribution to be 
homogenous and isotropic though time dependent). The velocity of photons through 
that medium is given by n = c (1 + (87rGp/3/c^) (1 — q/3)) Therefore if we assume 
that the photon velocity in an inertial frame does not exceed that flat space velocity 
c, then the equation of state of any form of matter p = qp must obey the strong 
energy condition q > 3. 


Higher Derivative Gravity: 


Drummond and Hathrell [3] obtained the higher derivative couplings which 
arise from the QED loop corrections to the graviton-photon vertex. The effective 
Lagrangian valid in the frequency range < a/ml is given by 


with the coefficients a = 


5 a ^ _ 26 o 


(45) 


^— 2 -,^ — STVT— 2 , c = (where a is the hne 

720 Trm^ ’ 620 Trm^ ’ 720 Trm^ ^ 

structure constant and rrif. the electron mass). The equations of motion from (45) are 
given by 


+ X^ = 0 


(46) 
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where 


= a + RV 

+2b - F-’'V^R\ - R\W^F^^^) (47) 

+4c + R^\r^f.F^n = 0 

Operating on equation (46) by V^ and using the Bianchi identity (10) the commutator 
identity (12) and the circular identity (13), we obtain the wave equation in higher 
derivative gravity: 


V^V/i-Fi/A + {Rp^uxF'^^ + R^~ 

+ [V,Xa - VaX,] = 0 


(48) 


In the derivations of [3-11] the terms in the curly bracket do not appear. This is be¬ 
cause in the derivations, the eikonal approximation is made in the Maxwell’s equation 
and therefore the Riemann and Ricci terms which arise on commuting the covariant 
derivatives do not appear. In (48) however we see that in the range of validity of the 
effective action (45), < a/ml, Riemann and Ricci terms in the curly brack¬ 

ets of (48) which are present in Einstein’s gravity are larger than the terms in the 
square brackets of (48) which arise from the loop corrections. The photon velocities 
are therefore given by the expression (31) for blackholes, and (44) for Friedmann- 
Robertson-Walker metric, and is subluminal even in the presence of higher derivative 
coupling terms. 
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